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A short introduction to the complex phenomena encountered in transition metal oxides with 
either charge or orbital or joint charge-and-orbital order, usually accompanied by magnetic order, is 
presented. It is argued that all the types of above ordered phases in these systems follow from strong 
Coulomb interactions as a result of certain compromise between competing instabilities towards 
various types of magnetic order and optimize the gain of kinetic energy in doped systems. This 
competition provides a natural explanation of the stripe order observed in doped cuprates, nickelates 
and manganites. In the undoped correlated insulators with orbital degrees of freedom the orbital 
order stabilizes particular types of anisotropic magnetic phases, and we contrast the case of decoupled 
spin and orbital degrees of freedom in the manganites with entangled spin-orbital states which decide 
about certain rather exotic phenomena observed in the perovskite vanadates at finite temperature. 
Examples of successful concepts in the theoretical approaches to these complex systems are given 
and some open problems of current interest are indicated. 
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I. DEGREES OF FREEDOM IN TRANSITION 
METAL OXIDES 

The physical properties of transition metal oxides are 
driven by strong electron interactions It is due to 
strong local Coulomb interactions that these systems ex- 
hibit very interesting and quite diverse instabilities to- 
wards ordered magnetic phases when doping x or tem- 
perature T is varied — is some cases also with orbital or- 
der. These instabilities are observed, inter alia, in rapid 
changes of the transport properties at the metal-insulator 
phase transitions, or in the onset of superconductivity. 

One of the outstanding problems in modern condensed 
matter theory is the description of strongly correlated 
electrons in various systems. When local Coulomb inter- 
actions are strong, the usual methods used for calculat- 
ing the electronic structure fail and have to be extended 
by the terms following from local interactions, either in 
the framework of the local density approximation (LDA) 
with Coulomb U, the so-called LDA+J7 method [2J, or 
by the self-energy within the dynamical mean-field the- 
ory (DMFT) |[, in the LDA+DMFT approach Q. This 
latter approach makes use of the local self-energy which 
becomes exact in the limit of infinite spatial dimension 
d = oo 5]. However, even these methods cannot over- 
come certain shortcomings of the effective one-particle 
theory which justifies modelling of these complex systems 
with Hamiltonians of the Hubbard type, and looking for 
solutions with methods of quantum many-body theory. 
The advantage of rapid progress in the electronic struc- 
ture calculations in recent years is that such models can 
nowadays use realistic parameters which follow from the 
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electronic structure calculations for a given system. 

Although the field of strongly correlated electronic sys- 
tems is very rich, we shall concentrate here on the phe- 
nomena observed in transition metal oxides. There are 
two major classes of systems with either perovskite struc- 
ture RMO3, or the layered structure R2MO4, with R 
standing for a rare-earth ion and M for a transition metal 
ion. In the latter class the subsequent layers of MOq oc- 
tahedra are displaced, so the electronic properties are 
well described by two-dimensional (2D) models, see Ref. 
1]. In both above structures electron correlations are 
strong and lead to remarkable consequences, with several 
degrees of freedom contributing simultaneously to coex- 
isting magnetic, charge and (in some cases also) orbital 
order. Examples of these complex phenomena are high- 
temperature superconductivity [qL the colossal magne- 
toresistance in the manganites |7 H9fl, a nd the Verwey 
transition in the magnetite (FesC^) |10j . Although the 
charge order occurs typically in doped systems, there are 
a few systems of formally mixed valence type, where the 
electron number per one transition metal ion is not an 
integer but local correlations stabilize charge order, as 
in the magnetite. The latter problem was recently ad- 
dressed and the mechanism of the Verwey transition was 
explained as triggered by the elcctron-phonon coupling 
enhanced by local Coulomb correlations [ll|, so we shall 
not discuss it here but refer an interested reader to an- 
other contribution in the same volume 

The electronic structure of transition metal oxides in- 
cludes several bands [l|, Ha ], but the properties of the 
system do depend on the states in the vicinity of the 
Fermi energy. It is usually sufficient to derive effective 
M-O-M hopping elements for a-bonds t a and 7r-bonds 
tit, and next use them in the effective model describing 
only 3d electrons [l4[. The respective kinetic energy is 
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described in a perovskite system by 
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Here {a, /?} = {x, z} are the indices of e g orbitals, 

\x) = (x 2 - y 2 )/V2, \z) = (3z 2 - r 2 )/V6 , (2) 

and this orbital flavor is in general not conserved along 
the hopping processes — the orbitals may be changed 
for the hopping along the bonds in ab planes in the per- 
ovskite structure. In contrast, the ti g orbital flavor, 

I a) = \yz), \b) = \zx), \c) = \xy). (3) 

is conserved for the hopping along the bonds in all three 
cubic directions 7 = a, b, c, as indicated by a single diag- 
onal hopping element t„ (for simplicity we assume only 
nearest neighbor hopping elements), with p, = {a, 6, c} 
labeling t 2g orbitals and referring to the cubic axes per- 
pendicular to the planes accomodating the respective or- 
bitals. The latter notation is introduced using an ideal 
cubic system in which each ti g orbital is perpendicular 
to a single cubic axis, for instance the \xy) orbital lies in 
the ab plane and is perpendicular to the c axis. 

On-site intraorbital Coulomb interactions are de- 
scribed by a single parameter U (identical for all 3d or- 
bitals): 
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In the simplest approach, the ratio U/W, where W is 
the bandwidth for the relevant partly filled band, decides 
whether electrons localize and the electronic structure 
changes to two Hubbard subbands in a Mott insulator, 
or the system is metallic, with rather strongly correlated 
electrons and possibly heavy effective masses (this hap- 
pens for the /-electron systems which are addressed in 
other contributions in this volume). As the hopping ele- 
ments along 7r bonds are significantly lower than the ones 
for g bonds [TiJ , the t 2g electrons in the early transition 
metal oxides (i.e. in titanium or vanadium oxides) are 
even stronger correlated than e g electrons in the -RM11O3 
or LaNiC>3 perovskites. This resembles the situation in 
molecular bonds in sp systems, with n bonds being al- 
ways stronger correlated than a bonds 15| . 

A second class of correlated insulators, so-called charge 
transfer insulators, arises when the oxyges states are 
within the gap between the two Hubbard subbands [l6| . 
A crucial parameter is the energy difference between the 
d and p electron (hole) levels, A = e p — £4 — here we 
use the hole notation relevant for the high-T c cuprates. 
When A > U one has a Mott-Hubbard insulator, but 
when A < U, the insulator is of charge transfer type. 

The electronic structure of the cuprates does not in- 
volve orbital degeneracy as the CuOg octahedra are elon- 
gated and the orbital degeneracy is removed for a tetrag- 
onal distortion. Therefore, s hole in the d 9 configuration 



occupies the |x) = (x 2 — y 2 )/^/2 orbital at each Cu 2+ 
ion in La2Cu04. The resulting charge transfer model for 
the Cu C>2 p lanes in the cuprates may be thus written as 
follows (l7| : 



Hd p = H + H n 



(5) 



H 



H.c. 



Hint — 



e p n p i tpd "/ m i (yd maa p lac 

tpp Vij (phaPjav + H - C -) . ( 6 ) 

(ij)aa 



The parameters of the charge transfer model ([5]) are: the 
oxygen energy s p (we assume that the reference d hole en- 
ergy Ed = 0), the d—p hybridization t p d, and the Coulomb 
interaction parameters for d and p orbitals, Ud and U p ; 
the same parameters describe also other Cu-0 s yste ms, 
as for instance CUO3 chains in YBa^Cu^Ofi-i- jr [181 ] . or 
CU2O5 coupled ladders in Sri4_ 2; Ca 3 ;Cu2404i [191 ]. Here 
n P i = n pl ^ + n P n and n pia = p\ a p la are charge density 
operators, j m i and r/ij are the phase factor for a pair of 
orbitals along the considered d — p (p — p) bond. The 
parameters for the cuprates which follow from the elec- 
tronic structure calculations are (in eV) 20]: A = 3.6, 
t pd = 1.3, tpp = 0.65, U d ^ 10.5, U d ^ 4.0. Electron cor- 
relations are moderate in spite of the large value of Ud 
[r^ ]. but they suffice to localize holes at Cu sites in the 
undoped system, such as La2Cu04 or YBa2Cu3C>6. Tak- 
ing the above parameters, A <§C Ud and these systems are 
charge transfer insulators, in contrast to the perovskite 
titanates and vanadates, which are Mott-Hubbard sys- 
tems. 

It is important to realize that the charge transfer gap A 
plays the role of an effective Coulomb parameter U = A 
in the correlated electronic structure of a charge transfer 
insulator. When the Cu-O-Cu hopping between two \x) 
orbitals along a bond in an ab plane is defined as t, this 
leads to the effective Hubbard model 21.1, 
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For the considered case of \x) orbitals the phase factors on 
each bond (ij) are identical. Note that t may be deduced 
from the charge transfer model ([5]), t = t 2 d / A = 0.4 eV, 
and for the actual ratio U/t = 10 the holes are strongly 
correlated. Hence, the undoped systems La2Cu04 or 
YBa2Cu30e are antiferromagnetic (AF) insulators. In 
general, the derivation of an effective model from the rel- 
evant multiband model is rather tedious — such a more 
complete model includes in addition next nearest (sec- 
ond) neighbor and third nearest neighbor hopping ele- 
ments {t',t"} and intersite Coulomb interactions [221 ]. 

A broad class of phenomena investigated for strongly 
correlated electron systems are the changes of their phys- 



3 



ical properties in the vicinity of metal-insulator transi- 
tions. As mentioned above, one way of localizing elec- 
trons in a correlated insulator is by changing the elec- 
tron interaction parameter U in Eq. © (or the charge 
transfer gap A). Although this may be easily realized 
only in theory, in certain systems the changes of the elec- 
tronic parameters are sufficient to induce metal-insulator 
transitions observed in V2O3 (23|. A more common sit- 
uation, however, is encountered in doped systems, where 
the carriers are released at certain doping concentration 
and the system becomes metallic. In contrast to the ear- 
lier suggestions, the one-band model is not sufficient to 
describe the metal- insulator transition in V2O3 (24[, and 
doping is not equivalent to varying external pressure [25[ ■ 
This and other metal-insulator transitions in the oxides 
are controlled by doping. A very well known example 
is the colossal magnetoresistance effect in the perovskite 
manganites [3] , another is the superconductivityin doped 
La2_ 2; Sr 2 ;Cu04 or YF^CuaOg+a; compounds [6]. Other 
examples can be found, for instance, in the excellent re- 
view article by Imada, Fujimori and Tokura [l[. 



In this paper we address in particular the phenom- 
ena related to magnetic and orbital order in transition 
metal oxides which follow from strong electron correla- 
tions. Charge order arises in doped systems, while the 
orbital order is common in transition metal oxides with 
partly filled degenerate orbitals. We begin in Sec. |TT]with 
the stripe phases in the cuprates, where we explain the 
stabilizing mechanism and show that the charge modula- 
tion is the way to optimize total energy in doped systems. 
While the properties of an undoped cuprate are driven by 
the AF superexchange, the systems with orbital superex- 
change interactions are more complex as the interactions 
are intrinsically frustrated [26J. These interactions are 
exemplified by the so-called compass model [13], see Sec. 
III[ and may give either highly degenerate ordered ground 
states, or the disordered orbital liquid. Consequences of 
the orbital superexchange for the magnetic order are ad- 
dressed in Sec. IIV| where we briefly summarize the struc- 
ture of the spin-orbital superexchange [28j . and demon- 
strate that spin and orbital degrees of freedom may be 
separated in the perovskite manganites. In contrast, in 
the perovskite vanadates spin-orbital entanglement plays 
a dominating role and decides about their properties at 
finite temperature, see Sec. |V] Finally, we give examples 
of coexisting charge-and-orbital order in doped systems 
in Sec. IVII A summary and some open problems in 
the field are given in Sec. I VIII Figures illustrating the 
theoretical concepts reviewed in this article will not be 
reproduced here — they may be found in the cited litera- 
ture which is far from being complete and was selected on 
the criterion of addressing the most important concepts 
in this field. 



II. STRIPE PHASES IN THE CUPRATES 

A crucial concept in the physics of the superconduct- 
ing cuprates is the Zhang- Rice singlet [2{|. It makes an 
explicit use of the charge transfer nature of the electronic 
structure, as a doped hole occupies not a Cu(d x ) orbital 
but a linear combination of p a orbitals with x 2 — y 2 sym- 
metry around a hole, which forms a singlet together with 
the hole at Cu ion. It is this concept which provides 
a justification for using the t—J model as the effective 
model describing the physical situation in the cuprates, 
and plays a prominent role in this class of compounds 

& 

The t-J model itself was derived from the Hubbard 



model in Cracow more than three decades ago 3lJ, us- 
ing the perturbation theory. A properly chosen canonical 
transformation leads from the full Hilbert space to an ef- 
fective low-energy Hamiltonian acting in the restricted 
space, where only spins and holes occur at different sites. 
It consists of the kinetic energy oc t and the superex- 
change interaction cx J between S — 1/2 spins: 

H t -j = -t (al*a ja +a] a h la )+J ^ fs r Sj + ^hifi^j , 
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with the superexchange interaction, 
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jected fermion operators and act in the restricted space. 
The above t—J model may also be derived directly from 
the charge transfer model — in this (realistic for the 
cuprates) case the superexchange includes both the An- 
derson and charge transfer excitations 32]. For the 
cuprates one finds J ~ 0.13 eV, which is either deduced 
from the magnetic experiments [33|, or derived from the 
charge transfer model using its parameters [2p] |. 

The first intriguing question concerning hole doping is 
whether a doped hole may propagate coherently in the 
antiferromagnet. Naively one might argue that a hole 
creates defects on its way, so it would need to make a 
hopping along a closed loop to annihilate these defects 
and to move in the square lattice with a minute disper- 
sion [34|]. Actually, this is the only process by which a 
hole may delocalize in the Ising model. The situation 
is quite different, however, when a hole is doped into a 
Heisenberg antiferromagnet — in this case the quantum 
fluctuations of the AF background may repair the defects 
created by the hole, and the hole dispersion occurs on the 
energy scale of J [35|, [36[ . This concept was confirmed 
by experiment, and indeed the hole dispersion on the low 
energy scale of J was observed in the cuprates [37j. De- 
tailed comparison between the experimental data of an- 
gle resolved photoemission experiments and the outcome 
of the theoretical calculations performed using the self- 
consistent Born approximation (SCBA) [36j were pre- 
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sented by several groups. Here we mention only the ex- 
perimental data of high quality obtained for S^CuC^C^ 
by Wells et al. [38| , which fit very well to the theoretical 
curves obtained with finite next neighbor hopping t' [39j j . 
This demonstrates that the right effective model for the 
high-T c cuprates is the t-t'-J rather than the t-J one. 

Higher doping of Cu0 2 planes leads to a gradual weak- 
ening of AF correlations, which however survive even in 
the overdoped regime at x ~ 0.2 (33|. There are sev- 
eral possibilities concerning the phase diagram of doped 
cuprates [9(, but the commonly accepted point of view 
now is that doped holes self-organize in form of phases 
with charge modulation [4(|. Such structures with co- 
existing charge and magnetic order, called stripe phases, 
were first discovered in the theory as an instability of 
doped antiferromagnets towards AF domains separated 
by (usually) nonmagnetic domain walls |4lj . Only a few 
years later their existence in the cuprates was confirmed 
in the neutron experiments of Tranquada et al. (42j. 
The stripe phases are characterized by the coexisting 
charge and magnetic order, with the charge density vary- 
ing twice faster than the spin density in the real space 

The first question concerning stripe phases is whether 
they would form as solitonic defects in the AF structure, 
i.e. in between different AF domains, or instead they 
are of polaronic nature not disturbing the AF order. Al- 
though naively one could argue that the polaronic mech- 
anism could give a better kinetic energy, this argument is 
misleading. To see this one can consider a cluster of three 
sites centered at the domain wall, filled by two electrons 
[43}. Due to strong correlations with U 3> £, the parti- 
cles are confined in this cluster, although it is just a part 
of the AF 2D plane. Taking two electrons with either 
identical spins or with opposite spins, it is straightfor- 
ward to estimate the ground state energy of two possible 
configurations: (i) polaronic (Ep) and (ii) solitonic (E$) 
one. One finds that the solitonic energy is lower by the 
superexchange energy J which arises from the three-site 
hopping terms in this cluster (43[, with: Ep = —\f2t 
and Eg = — y/2t — At 2 /U. This simple argument explains 
the experimental finding that charge walls separate AF 
domains with different phase of the order parameter. 

Quantitative results for the stripe phases were first 
obtained using the Hartree-Fock (HF) approximation 
[4ll |4"3| . and then refined using variational wave func- 
tions Hi!], within DMFT for the stripe phases 0, HtJ, 
and slave-boson approach [48[ . As usually, the HF serves 
only as a hint for possible instabilities, and gives remark- 
ably robust stripe structures pij with rather large am- 
plitude of the charge density between the domain wall 
and the centers of AF domains, and the filling of half of 
doped hole per one stripe charge unit cell, as observed 
in experiment. These stripes are vertical (or horizontal), 
meaning that domain walls are along (10) [or (01)] di- 
rection, and insulating. Actually, their stability follows 
from a small gap which opens in the electronic structure. 
This mechanism is subtle and involves certain additional 



modulation, either spin or charge density wave, along the 
domain walls [3], so one has to expect major changes 
when electron correlations are implemented. However, 
variational calculations confirmed this picture to some 
extent [45|j , although the question whether the stripes 
are insulating or not could not be resolved (following the 
HF results, it was believed for a long time that the stripes 
are insulating). 

Stable stripes were also found using an exact di- 
agonalization method within the DMFT for the two- 
dimensional Hubbard model, in the broad doping range 
0.03 < x < 0.2 in La 2 _ a; Sr a; Cu04 These calculations 
allowed also to reproduce the observed crossover from di- 
agonal (ll]to vertical (01) site-centered stripes at doping 
x ~ 0.05 j49|. In addition, also the doping dependence of 
the size of magnetic domains and chemical potential shift 



A/i cx 



were found to be in quantitative agreement 



with the experimental results for L^-^Sr^CuO.^ In this 
way the paradigm of insulating stripe phases was abol- 
ished — the chemical potential was varying with doping 
within the metallic phase. 

The spectral functions obtained within the DMFT (47j 
show a coexistence of the incoherent states in the lower 
Hubbard band and a coherent quasiparticle (QP) close to 
the Fermi energy. The main features of the spectra are: 
a flat part of the QP band near the X ~ (tt,0) point, 
and gaps for charge excitations at the Y = (0, ir) and 
S = (it/2, 7r/2) points in the low-doping regime x < 1/8. 
These gaps are gradually filled and close under increas- 
ing doping, in agreement with the experimental data for 
La2- a; Sr a ;Cu04 obtained using angle resolved photoemis- 
sion [50j . In the range of low temperatures T the obtained 
spectra have a distinct QP peak at the X — (n, 0) point, 
present just below the Fermi energy fj,, and a charge gap 
and well defined QP at the S = (tt/2,tt/2) point [471 ] . 
These calculations demonstrated the importance of dy- 
namical correlations which strongly screen the local po- 
tentials resulting from on-site Coulomb interactions and 
lead thus to drastic changes in the distribution of spec- 
tral weight with respect to the HF picture. It was also 
shown that the melting of stripe order is influenced by the 
second neighbor hopping element t' , which plays also an 
important role in explaining the observed difference in 
the spectral properties between Bi 2 Sr 2 CaCu 2 08+y (5l| 
and La2- a: Sr a: Cu04 [5(j. At the same time, t' can tip 
the energy balance between the filled diagonal and half- 
filled vertical stripes [481 ] . which might explain a change 
in the spatial orientation of stripes observed in the high 
T c cuprates at the doping x ~ 1/16. 

More insight into the charge and magnetization distri- 
bution as well as into the stability of stripe phases could 
be obtained using a rotationally invariant version of the 
slave-boson approach in spin space in the 2D Hubbard 
model [48J. This approach allowed one to treat strong 
electron correlations in the stripe phases with large unit 
cells relevant in the low doping regime, and gave results 
representative of the thermodynamic limit. It also helped 
to resolve the longstanding controversy concerning the 



5 



role played by the kinetic energy in stripe phases. While 
the transverse hopping across the domain walls yields 
the largest kinetic energy gain in the case of the insu- 
lating stripes with one hole per site, the holes propagat- 
ing along the domain walls stabilize the metallic vertical 
(01) stripes with one hole per two sites, as found in the 
cuprates. 

Recently observed pattern of unidirectional domains 
in high-T c superconductors [52| motivated also search 
for coexisting charge modulation and d-wave supercon- 
ductivity. Indeed, half-filled charge domains separated 
by four lattice spacings were obtained along one of the 
crystal axes leading to modulated superconductivity with 
out-of- pha se d-wave order parameters in neighboring do- 
mains (53|. Both renormalized mean- field (MF) theory 
and variational Monte Carlo calculations yield that the 
energies of modulated and uniform phases are very close 
to each other, so modulated phases could easily be sta- 
bilized by other effects going beyond the t—J model used 
in these calculations. Novel doped phases with supercon- 
ductivity coexisting with charge modulation or even the 
valence-bond solid order were also reported recently . 



III. INTRINSIC FRUSTRATION OF THE 
ORBITAL SUPEREXCHANGE 

In the undoped transition metal compounds the 
physics is however frequently not so simple as in CuC-2 
planes of high T c superconductors, where the superex- 
change stabilizes the AF long-range order. This happens 
in particular when the orbital degrees of freedom are ac- 
tive (when degenerate orbitals are only partly filled) and 

. The central 
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contribute to the magnetic order 
property of the orbital degrees of freedom is that they 
are intrinsically frustrated [HI, H?} , so they may lead to 
novel (ordered or disordered) phases. 

Frustration in magnetic systems may be of geometrical 
origin if only nearest neighbor interactions are present, 
or may arise due to competing exchange interactions 
[HI, [59[ ■ For instance, when one considers FM interac- 
tions along every second vertical line in the square lattice 
while all other interactions are AF — then this 2D Ising 
model is exactly solvable and has a lower critical tem- 
perature [60l ] than the one with isotropic exchange inter- 
actions. Frustration for quantum spins acts to enhance 
the effects of quantum fluctuations, leading to a num- 
ber of different types of magnetically disordered states, 
among which some of the more familiar ones are static 
and resonating valence-bond (VB) phases (59|. However, 
also ordered phases may emerge in systems with frus- 
trated spin interactions from their disordered manifolds 
of states, and their mechanism of stability is nowadays 
called "order-by-disorder" [HH . Numerous materials are 
known at present whose physical properties could be 
understood only by employing microscopic models with 
frustrated spin interactions in which some of these theo- 
retical concepts are exemplified. 



A prototype model to study frustration in pseudospin 
systems which mimic the directional orbital superex- 
change [6l| is the 2D compass model (27| 
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In this model the rf rj interactions oc J x for horizontal 
bonds (ij) (along the a axis) compete with the t?t? ones 
cx J z in the vertical direction (along the b axis). Re- 
cently the structure of eigenstates in this model was in- 
vestigated by numerical methods |62j ]. and it was shown 
using quantum Monte Carlo that a phase transition at 
finite temperature exists in the 2D compass model [63| . 
suggesting that this model is indeed in the 2D Ising uni- 
versality class. A competition of pseudospin interactions 
along different directions results here in intersite cor- 
relations similar to the anisotropic XY model, and in 
competition between two types of Ising-like order. This 
competition culminates in the highly degenerate ground 
state at the compass point (i.e., when all interactions 
have the same strength) [62|, and generates there a first 
order phase transition when the anisotropic model with 
Jz > Jx changes into J z < J x through the J z = J x 
transition point [64j . It is interesting to note that a sim- 
ilar first order quantum phase transition occurs also in 
the one-dimensional (ID) compass model [65j . when both 
above interactions alternate along the chain (N' = N/2 
is the number of unit cells): 
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n 1D = E{ J -' 
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This model was solved exactly in the entire range of 
{Jx,Jz} parameters [65| by mapping onto the exactly 
solvable quantum Ising model |66j in different subspaces. 
Equal coupling constants J x = J z = J correspond here 
to the the quantum critical point, where the disordered 
phase (orbital liquid) emerges from two different types of 
hidden order, and the first order transition takes place. A 
similar transition was shown as well by an exact solution 
of the compass ladder ■ 

The compass model is currently under discussion also 
due to its interdisciplinary character. It can be derived 
using the symmetry arguments which are necessary for 
the realization of doubly degenerate states which are pro- 
tected from external perturbations in a wide class of 
Hamiltonians [68j . The Hamiltonian with this symmetry 
can be physically implemented in Josephson junctions, 
and it was argued that these junctions provide fault tol- 
erant quantum bits. Recently magnetic interactions in 
Mott insulators with strong spin-orbit coupling were also 
discussed and it was pointed out [iU that they may pro- 
vide a realization of the exactly solvable Kitaev model 
on the honeycomb lattice, which is relevant for quantum 
computation [70| . 



6 



IV. SPIN-ORBITAL SUPEREXCHANGE 

In several transition metal oxides with active degrees 
of freedom one finds coexisting magnetic and orbital or- 
der, both in Mott and in charge transfer insulators. Ex- 
perimental observations give frequently some unexpected 
properties and the question arises how to explain partic- 
ular types of observed coexisting spin-orbital order. First 
of all, magnetic interactions in these systems frequently 
break the cubic symmetry of the perovskite lattice and 
AF phases arise with ferromagnetic (FM) interactions 
along certain directions (in contrast to the G-AF phase 
with isotropic AF interactions, the same along each cu- 
bic direction). Two of them are quite common: (i) the 
C-AF phase with FM interactions along the c axis as in 
LaVC>3, and (ii) the A-AF phase with FM interactions 
within the ab planes as in LaMn03 [57[. These phases 
follow from the microscopic models (see below) which jus- 
tify the complementary behavior of the observed (orbital 
and magnetic) order postulated by Goodenough in the 
manganites [71|: alternating orbital (AO) order supports 
FM spin order, while ferro-orbital (FO) order supports 
AF spin order. The structure of spin-orbital superex- 
change described below allows to understand better the 
physical mechanism beyond this complementarity, known 
since long as the Goodenough-Kanamori rules [72| . 

Realistic superexchange models for transition metal 
oxides with orbital degrees of freedom contain both spin 
and orbital operators. They may be derived by consid- 
ering intersite charge excitations in a Mott insulator in 
a way similar to the derivation of the t—J model from 
the Hubbard model (and applying to the cuprates), as 
described in Sec. [TTJ This task is somewhat involved, so 
usually one considers an effective model with hopping ele- 
ments between d orbitals of transition metal ions derived 
from electron transitions over the intermediate oxygen 
orbitals [l4| , as in Eq. ((T|) . The energy scale for the hop- 
ping is set by the largest hopping element t: the (dda) 
element in case of e g systems, and the (ddn) element 
when only 7r electrons contribute in systems with degen- 
erate and partly filled <2g orbitals. For noninteracting 
electrons the Hamiltonian Hq (with crystal-field terms) 
would lead to tight-binding bands, but in a Mott insu- 
lator one is in the regime of large Coulomb interaction 
U >■ t, so charge fluctuations are suppressed and the 
hopping elements can only contribute via virtual excita- 
tions, leading to the superexchange as described below. 

The derivation of the superexchange involves virtual 
charge excitations d^d" 1 ^ d™ 4 " 1 ^ -1 between two 
neighboring transition metal ions with m 3d electrons 
each, which have to be included with the correct ex- 
citation energies following from the structure of local 
Coulomb interactions in degenerate d states (a similar 
expression can also be written for / electrons) . These in- 
teractions are well known since long [73| , but nevertheless 
simplified expressions can be still found in the literature 
which in some cases lead to misleading or even qualita- 
tively incorrect conclusions. When only one symmetry 



class of electrons is involved in charge excitations, either 
e g or t2 g electrons in a perovskite system, the interactions 
read: 



Hint — U 



^2 n iaf n ial + (U' - - Jh) ^ 



i,a</3 

2Jh ^ • Sip , 

i,a<(3 



(13) 



with a = —a. The parameters {U, Jh} (with U' = 
U — 2Jh) determine the excitation energies by the cor- 
responding eigenenergies of the local electron-electron 
Coulomb interaction Hamiltonian (fl~3|) . The eigenstates 
have to be derived for each particular situation sepa- 
rately. For instance, if two-hole d 8 excited states are con- 
sidered in the d^d^ ^ d s dj° transitions for the KCuF 3 
cuprate, there are triplet (S = 1) and singlet (S = 0) 
states when two holes in the d s configuration occupy dif- 
ferent orbitals, with the energies U — 3Jh and U — Jh, 
as well as two other (intraorbital) singlet states with the 
energies U — Jh and U + Jh [ZJ] • Note that a double oc- 
cupancy of either e g orbital (with energy U), considered 
in the context of the intraorbital Coulomb interaction 
P5p is not an eigenstate of the local Hamiltonian (fT3"|) . 

In a general case the interorbital interactions, Coulomb 
U a /3 and exchange J a p elements replace U' and Jh — in 
contrast to the intraorbital ones (|15[) they are anisotropic, 
but satisfy a constraint which guarantees the invariance 
of interactions in the orbital space [75j , 



U — U a p + 2J a p 



(14) 



for each pair of interacting orbitals {a/3}. These interac- 
tions are frequently parametrized by the Racah parame- 
ters {A,B,C} [76|, and one finds that 



U = A + 4B + 3C , 



(15) 



and is identical for all 3c? orbitals, while Hund's exchange 
depends on the orbital states of the pair of interacting 
electrons 



for e g electrons, 
J H =4B 



C. 



and it is somewhat stronger than for ti g electrons, 
J H = 3B + C . 



(16) 



(17) 



Thus the parameter Jh used in Eq. (flU]) refers to the 
above values depending on whether a system with e g or 
ti g orbital degrees of freedom is considered. More details 
about the structure of local Coulomb interactions which 
depends on a single parameter, 



>1 : 



Jh 
U 



(18) 



may be found in Ref. [75| , while the experimental values 
of the Racah parameters {B,C} which are known with 
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high accuracy from the atomic spectra are given in Ref. 
[77j for several transition metal ions. Unfortunately, the 
value of U (or A) is known only with much lower accuracy 
and the mechanism of screening which leads to the values 
U ~ 5 — 10 eV is difficult to implement in the theory. 
Hence, the value of U is frequently used as a parameter, 
unless it can be derived from the experimental data, as 
for instance from the optical excitations, see below. 

The above structure of the electron-electron interac- 
tions (|13[) determines the excitation energies e n due 
to the multiplct states in charge transitions d™d™ f± 

dj 1 which enter the superexchange in the respec- 
tive denominators of 4t 2 /e n . Examples of such spectra 
are presented in Ref. [28| • As a rule, the high-spin states 
have the lowest energy U—3Jh independently of the elec- 
tron number to in the electronic configuration d m under 
consideration, while the energies of low-spin states de- 
pend on to, and may even contain fractions of Jjj due 
to the anisotropy of J a p Hund's elements, as found for 
instance in the case of LaMnC>3 [781 ]. 

The spin-orbital superexchange is the effective low- 
energy Hamiltonian which involves products of spin and 
orbital operators. The spin interactions are described by 
spin scalar products Si ■ Sj on each bond (ij) connect- 
ing two nearest neighbor transition metal ions and obey 
the SU(2) symmetry, while the orbital operators {f,fj}, 
with fi — {t^ , Tj~ , r*}, obey only much lower symmetry 
(at most cubic for a cubic lattice) and appear either as a 
scalar product f • fj, or only as certain components, for 
instance an Ising term rfrf, resembling the interactions 
in the 2D compass model (flTj) . In general the spin-orbital 
superexchange is of the form (28j . 



n., 



J H {j^fi-Si + s^+icV}, (i9) 



<y>ll7 



with the constant J defined in Eq. (fT0| . The orbital op- 
erators and it^ depend on the direction 7 = a, b, c 
in the cubic lattice and involve the active orbitals on each 
bond (ij) (either e g or t2 g ) along direction 7 — they par- 
ticipate in d^d" 1 ^ <i™ +1 d™ _1 virtual excitations, and 
thus these interactions have the symmetry of the lattice 
(i.e. cubic symmetry in the perovskites). As an example 
we introduce here the superexchange between V 3+ ions 
in the d 2 configuration with S = 1 spins p79| , as realized 
in .RVO3 perovskites considered below in Sec. [V] and 
controlled by the orbital operators: 



1 



rjr 3 Ti x tj + -n^n 



-nr^m + rij)} (20) 



(7) 



4 / 2 



k[f = { rjrx [ fi ■ fj + -n l n j ) + r]r 3 [ f { x f } + -n^j 



j(l + Vn)(ni +rij) 



(7) 



(21) 



They arise from the d 2 d 2 ^ dfdj charge excitations, 
leading either to high-spin or to low-spin configura- 
tions, so Hund's exchange in the multiplet structure of 
a V 2+ ions enters via the coefficients r*i = 1/(1 — 3?y) 
and r 3 = 1/(1 + 2ry) (the low-spin excitations occur also 
at energy U, so the corresponding coefficient is T2 = 1). 
Eqs. (|2^|) and (|2T|) are general and refer to two active 
orbital flavors along the cubic axis 7. The leading or- 
bital interactions are proportional to the scalar products 
(fi ■ fj)^ of orbital operators on the bonds as both or- 
bitals are active and may generate charge excitations, 
but the structure of local Coulomb interactions (p~3|) is 
responsible for additional terms, 



(fi X fj) 



(c) 



{-rfrf +T i T j )+- 



(22) 



which violate the conservation of the orbital quantum 
numbers. The operator stands for the number of 
active electrons at site i along the bond (ij), for instance 

...(c) 

for a bond along the c axis this number is n\ — n ia + n ib 
[this notation for the ti g orbitals is defined in Eq. 

The superexchange model (TT91 consists typically of 
several terms which originate from different charge exci- 
tations. This feature made it possible to relate the aver- 
ages of these different excitations to the spectral weights 
in the optical spectroscopy [80j , and serves now as a stan- 
dard theoretical tool to explain the observed anisotropy 
and temperature de pen dence of the spectral weights in 
the optical spectra 28]. In a correlated insulator the 
electrons are almost localized and the only kinetic en- 
ergy which is left [8l| is associated with the same virtual 
charge excitations that contribute also to the superex- 
change (|19[) . Therefore, we one may define the individual 
kinetic energy contributions Kn\ which refer to differ- 
ent energy regimes in the optical transitions and can be 
determined from the superexchange using the Hellman- 
Feynman theorem [82| , 



KM = -2(H^(ij)). 



(23) 



(7) 



For convenience, we define the energy contribution K 
for the Hubbard subband n as a positive quantity. 

The magnetic properties of the transition metal oxides 
with active orbital degrees of freedom are usually dis- 
cussed in terms of magnetic exchange constants which de- 
termine both the type of the magnetic order in the ground 
state (at T = 0) and the magnetic excitations (magnons 
observed in the neutron scattering experiments). The ex- 
change constants are usually found for a bond (ij) along 
each nonequivalent axis 7 by averaging over the orbital 
operators in Eq. (fl9|) . 



Jij — (3 i 



ij I 



(24) 



which leads to an anisotropic spin exchange model Hamil- 
tonian adequate for anisotropic magnetic phases, such as 
for instance A-AF or C-AF phase realized for instance in 
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LaMn0 3 and LaV0 3 , 



H s — J a 



Si • Sj -\- J c 



Si ■ Sj 



(25) 



This procedure assumes implicitly that spin and or- 
bital operators can be separated from each other and 
ignores the possibility of quantum entanglement [83| and 
of c omp osite spin-orbital excitations introduced in Refs. 
[74l |84| . It turns out that such excitations play a promi- 
nent role in destabilizing the classical AF long-range or- 
der in the d 9 spin-orbital model |26j , and have observable 
consequences in the perovskite vanadates, see Sec. |V] 

In some cases, however, the spin and orbital degrees of 
freedom may be disentangled and the obtained theoreti- 
cal results explain well the experimental findings. One of 
the best examples is LaMnOs, where the exchange con- 
stants {J a b, Jc} deduced from the neutron scattering (8f| 
can be explained by the superexchange model assuming 
a classical ansatz for the ground state with AO order, 



i*o> = n 



i£A 



n 



■>B)j, 



(26) 



with the orbital states, \0A)i and |6*s)j, characterized by 
opposite angles (9 a = —9 b) on two sublattices A and B 
in the ab planes, and repeated in the subsequent planes 
along the c axis. The AO order is stable below the orbital 
transition temperature Too — 760 K [90| , which is rather 
high compared with the Neel temperature T/v — 140 K 
— therefore one may consider the AO order between the 
sublattices (i € A, j € B), 



\0 



A/i 



cosl- 



4 + 



am- 



i- 



cus 



2 



•ami- 



13 ■ 



(27) 



as frozen in the temperature range T < 300 K rele- 
vant for the magnetic excitations. Here we used the 
simplified notation for the e g orbital basis introduced 
in Eq. (j2J|. Using the well motivated parameter set, 
the experimental values of the exchange constants in 
LaMnOa are reproduced by the angle 9 ~ 94°, but 
it has been shown that somewhat higher values of 9 
would also be consistent with a model including explic- 
itly charge transfer superexchange terms [28| . In any 
case, the occupied orbitals in LaMnOs are closer to sym- 
metric/antisymmetric combinations of {\z}, \x)} than to 
the directional 3x 2 — r 2 /3y 2 — r 2 orbitals, as illustrated 
in the early literature on the subject. 

The experimental proof that the spin and orbital op- 
erators may be disentangled in LaMn03 is provided by 
the optical spectroscopy, which shows rather distinct 
anisotropy of the low-energy spectral weights between the 
polarization in ab planes on the one hand and along the 
c axis on the other (86|. It is quite remarkable that the 
temperature dependence of these spectral weights may 
be very well explained by the spin-orbital superexchange 



model of Ref. [78j, using the same parameters as those 
used to calculate the exchange constants, in the broad 
temperature range < T < 300 K [86| . It follows alone 
from the temperature variation of spin correlation func- 
tions for a bond (ij) within the ab planes and along the 
c axis, 



Sab = (Si ■ Sj)^ 



= (^-^> (C) 



(28) 



while the AO order remains unchanged, as explained in 
Ref. 

In spite of the disentangled spin and orbital dynamics 
in the i?MnOs perovskites, where i?,=Lu,- • - ,La stands for 
a rare earth atom, their several properties are not fully 
understood. One of them is the nature of the insula ting 
state which comes partly due to JT interactions |87H89j 
and partly due to the orbital superexchange interactions 
[7^ | . Another puzzling feature is the phase diagram of the 
i?Mn03 family of compounds, where i?=Lu,- • -,La stands 
for a rare earth atom, which exhibits a phase transition 
from the A- AF to a rather peculiar .E-AF phase [9(| ■ 



V. FINGERPRINTS OF SPIN-ORBITAL 
ENTANGLEMENT IN THE i?VO a PEROVSKITES 

A. Spin-orbital entanglement 

The coupling between spin and orbital operators in the 
spin-orbital superexchange may be quite strong in some 
cases — the excellent example of this coupling are the 
vanadium perovskites, see below. Although the C-AF 
phase observed in the entire family of -RVO3 compounds 
[9l|,[92|], where i?=Lu,- • -,La stands for a rare earth atom, 
satisfies to some extent the Goodenough-Kanamori rules 
[72|, with FM order along the c axis where the active a 
and b orbitals ^ alternate — the AO order is very weak 
here and the orbital fluctuations play a very important 
role [Z9|]. This situation is opposite to the frozen AO 
order in LaMnOs, which can explain both the observed 
magnetic exchange constants and the distribution of the 
optical spectral weights. In LaV03 the FM exchange in- 
teraction is enhanced far beyond the usual mechanism 
following from the splitting between the high-spin and 
low-spin states due to finite Hund's exchange Jh- Evi- 
dence of orbital fluctuations in the i?V03 perovskites was 
also found in pressure experiments, which show a distinct 
competition between the C-AF and G-AF spin order, ac- 
companied by the complementary G-AO and G-AO order 
of {a,b} orbitals pf . 

To understand better the essence of entangled spin- 
orbital states, we present first the results of the model 
calculation with four-site chains along the c axis, de- 
scribed by the spin-orbital superexchange models rele- 
vant for titanates and vanadates. These calculations 
served to identify spin-orbital entangled states for in- 
creasing multiplet splitting oc 77 (83|. A prototype model 
to study frustration and entanglement in coupled spin 



9 



and pseudospin (orbital) systems is the one-dimensional 
(ID) SU(4) model [93|. This example is remarkable, as in 
a purely spin ID model one expects no frustration when 
only nearest-neighbor interactions are present. However, 
both spins and pseudospins appear here on a completely 
symmetrical and equal footing with joint spin-pseudospin 
operators and compete with each other, forming a group 
of elementary generators in the SU(4) symmetry. Three 
types of elementary excitations contribute to the thermo- 
dynamic properties: spin, orbital, and joint spin-orbital 
ones. This is indeed confirmed by the entropy data of this 
model obtained from a numerical analysis 94] , which in- 
creases three times faster than that of the ID AF Heisen- 
berg model. This also implies that the intersite correla- 
tion functions are intimately interrelated and it is impos- 
sible to separate the two subsystems. Hence, one has to 
treat explicitly entangled spin-pseudospin states. 

As a useful tool to verify the Goodenough-Kanamori 
rules [72| in spin-orbital models for 7j 2g electrons with 
spins either S =1 or S =1/2 and pseudospins r = 1/2, 
we introduce spin and orbital correlations defined for a 
bond (ij), 

S i j = {8 i -Sj)/(2S)\ T ij = {f i -f j ). (29) 

When they are compared with each other, and with the 
composite spin-orbital correlation function defined as a 
difference between the exact value and the MF factorized 
correlations on a bond (ij) [8a |. 

Cij = {((Si ■ ■ fj)) - (5, • SjXfi ■ fj)} /(2S) 2 , 

(30) 

one may conclude whether the spin and orbital operators 
are disentangled. If Cij = 0, the spin and orbital opera- 
tors are disentangled and their MF decoupling is exact, 
while if Cij < — spin and orbital operators are entan- 
gled, and joint spin-orbital fluctuations contribute even 
at T = 0. Two spin-orbital models were investigated in 
Ref. [HI : (i) the titanate model for d 1 ionic configura- 
tions of Ti 3+ ions in the i?TiC>3 perovskites with S =1/2 
[95j . and (ii) the vanadate model for d 2 configurations of 
V 3+ ions in the RV0 3 perovskites with S = 1 0. For 
more details about the structure of the superexchange 
7~Lj ([H]) in both models see for instance Ref. 28]. 

As the chain-like cluster is ID and only two orbital 
{a, b} flavors contribute in each case, one recovers the 
SU(4) model in the d 1 (titanate) case at 77 = 0, and 
Sij = Tij = Cij = -0.25 for TV = 4 sites % a 

closer inspection one finds that the ground state wave 
function for the four-site cluster is close to a total spin- 
orbital singlet, involving a linear combination of (spin 
singlet/orbital triplet) and (spin triplet /orbital singlet) 
states for each bond (ij). This result manifestly contra- 
dicts the celebrated Goodenough-Kanamori rules [72[ , as 
both spin and orbital correlations have the same sign. 
At finite rj the SU(4) degeneracy of all intersite corre- 
lations is removed — one finds T,-j < C,-j < < 
in the regime of spin singlet (S — 0) ground state, 



and the Goodenough-Kanamori rule with complementary 
spin/orbital correlations is still violated. A qualitatively 
similar case is found in a mathematical SU(2)(g>SU(2) 
model (not realized in transition metal oxides) [96] , where 
the ground state is entangled in a broad range of param- 
eters, including the exactly solvable case with alternating 
spin and orbital singlets on the bonds }97j . 

The vanadate d 2 model behaves also in a similar way in 
a range of small values of 77, with all three Sij , Tij and Cy 
correlations being negative. Again, the composite spin- 
orbital correlations are here finite (Cy < 0), spin and or- 
bital variables are entangled, and the MF factorization of 
the ground state into spin and orbital part fails. Only for 
sufficiently large r\ do the spins reorient in the FM ground 
state, and decouple from the orbitals. In this regime, cor- 
responding to the experimentally observed C-AF phase 
of LaV03 (and other cubic vanadates), spin-orbital en- 
tanglement ceases to exist in the ground state. However, 
as we will see below, it has still remarkable consequences 
at finite temperature, where entangled spin-orbital states 
again play a role. 

A crucial observation concerning the applicability of 
the Goodenough-Kanamori rules to the quantum models 
of t2 g electrons in one dimension can be made by com- 
paring spin exchange constants Jy calculated from Eq. 
(f2"4")l with the actual values of intersite spin correlations 
Sij (12"9"|) . One finds that exchange interaction is negative 
(Jij < 0), so formally favors FM spin orientation, in the 
singlet phase at low values of rj, but it is accompanied 
by AF spin correlations (Sij < 0). This demonstrates 
that the ground state energy calculated in the MF the- 
ory would be enhanced, so the MF approach cannot be 
used [83[ . This result follows from large spin-orbital fluc- 
tuations which cause also large fluctuations of the ex- 
change constants around the average value, measured by 

SJ = {((Ji^f) ~ JfjY /2 - Altogether, this result chal- 
lenges the usual interpretation of the magnetic data in 
the spin-orbital systems with the exchange constants de- 
termined by averaging over the orbital operators, see Eq. 
(f2"4")l . Fortunately, Hund's exchange is large enough in 
real materials and this conceptual difficulty is removed 
in transition metal oxides, but one may expect that ex- 
perimental results in the range of finite temperature will 
depend on the entangled states discussed above. 



B. Phase diagram of the KVO3 perovskites 

The phase diagram of the i?V03 perovskites [9l], [92j 
is qualitatively different from the one for the -RVO3 per- 
ovskites [90( | and indicates the proximity of spin and or- 
bital energy scales. Experimental studies have shown 
that the C-AF order is common to the entire family of 
the -RVO3 vanadates, and in general the magnetic tran- 
sition occurs below the orbital transition, ?Vi < Too, 
except for LaV0 3 with T N1 ~ T 00 EHhI- When the 
ionic radius rn decreases, the Neel temperature Tjvi also 
decreases, while the orbital transition temperature Too 
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increases, passes through a maximum close to YVO3, and 
next decreases towards LUVO3. This provided an experi- 
mental challenge to the theory which was addressed only 
recently using the spin-orbital superexchange model [98| . 
One hnds that the C-AF order develops in LaVC>3 be- 
low T/vi ~ 143 K, and is almost immediately followed by 
a weak structural transition stabilizing the weak G-AO 
order at T o — 141 K [H HI]. This provides a con- 
straint on the theoretical model. Remarkably, the mag- 
netic order parameter in the C-AF phase of LaVC>3 is 
strongly reduced to ~ 1.3/is, much below the reduction 
expected from quantum fluctuations in the C-AF phase 
(being only 6% for S = 1 spins [§9j]) — also this reduction 
of the measured magnetization could not be explained so 
far. 

In order to unravel the physical mechanism responsi- 
ble for the decrease of Too from YVO3 to LUVO3 one 
has to analyze in more detail the evolution of GdFeC>3 
distortions with for decreasing ionic radius rn 98l| . Such 
distortions are common for the perovskites [lOOj , and one 
expects that they should increase when the ionic radius 
tr decreases, as observed in the i?MnC>3 perovskites [90| . 
In the -RVO3 family the distortions are described by two 
subsequent rotations of VOq octahedra: (i) by an angle 
$ around the b axis, and (ii) by an angle ip around the 
c axis. Increasing angle d causes a decrease of V-O-V 
bond angle along the c direction, being ir — 2'd, and leads 
to an orthorhombic lattice distortion u — (b—a)/a, where 
a and b are the lattice parameters of the Pbnm structure 
of -RVO3. By the analysis o f the structural data for the 
i?V0 3 perovskites [Toil. [l02 1 one finds the following em- 
pirical relation between the ionic radius tr and the angle 
§: 



tr = r - a sin 1? , 



(31) 



where ro = 1.5 A and a = 0.95 A are the empiri- 
cal parameters. This allows one to use the angle $ to 
parametrize the dependence of the microscopic parame- 
ters of the Hamiltonian and to investigate the transition 
temperatures Too and Tjvi as functions of tr. 

The spin-orbital model introduced in Ref. [98[ to de- 
scribe the phase diagram of -RVO3 reads: 



{ij)\\c (ij)\\ab 

;Y,T* + \NK{u-u Q (d))\ 



'iQ T z 



gu 



(32) 



where 7 = a, b, c labels the cubic axes, and the opera- 
tors are given by Eqs. (|2"0")) and (|2"Tj) . The superexchange 
is supplemented by the crystal field term cx E z , the or- 
bital interactions terms oc V c and cx V a b induced by lat- 
tice distortions, and the orbital-lattice term oc g which 
is counteracted by the lattice elastic energy oc K. All 



these terms are necessary in a realistic model which re- 
produces the behavior of the -RVO3 perovskites at finite 
temperature. 

The crystal field splitting breaks the cubic symmetry 
in distorted VOg octah edra , as obtained in the electronic 
structure calculations 103] and from the point charge 



model [98|, and the actual filling of t 2g orbitals is: 



= 1, 



riia + n ib = 1 , 



(33) 



so the superexchange (|2U)) and (|2"Tj) in Eq. (|3"2"|) is ex- 
pressed by the orbital operators 73 = {t,- ,t^~ , rf} (and 
their components) as explained in Sec. IV. The splitting 
oc E z between a (yz) and b (zx) orbitals is given by the 
pseudospin rf operators, 



(34) 



which refer to two active orbital flavors {a, 6} in RVO3. It 
is characterized by the vector Q = (tt, n, 0) in reciprocal 
space — it alternates in the ab planes, but is uniform 
along the c axis. Thus, this splitting competes with the 
(weak) G-AO order supporting the observed C-AF phase 
at temperature T < Tjyi. 

In addition, the model (I3"2l includes: (i) intersite or- 
bital interactions oc V a b, V c (which originate from the cou- 
pling to the lattice), and (ii) orbital-lattice term oc g 
which induces orbital polarization (rf) 7^ when the 
lattice distortion u increases. The orbital interactions 
induced by the distortions of the VOg octahedra and by 
GdFe03 distortions of the lattice, V a b > and V c > 0, 
also favor the C-AO order (like the crystal field E z > 0). 
Note that V c > counteracts the orbital interactions in- 
eluded in the superexchange via K\y operators ([21]). The 
last two terms in Eq. (1321) describe the linear coupling 
oc g > between active {yz, zx} orbitals and the or- 
thorhombic lattice distortion u. The elastic energy which 
counteracts lattice distortion u is given the force constant 
K, and N is the number of V 3+ ions. The coupling oc gu 
acts as a transverse field in the pseudospin space. While 
the eigenstates ^(|a) ± favored by rf cannot be re- 
alized due to the competition with all the other terms, 
increasing lattice distortion u (increasing angle $) modi- 
fies the orbital order and intersite orbital correlations. 

The crystal field splitting E z ($), orbital interactions 
{Vab{fl), Vcft)}, and the orbital-lattice coupling g e ff(i?) = 
gu depend on the tilting angle 1). In case of V c one 
may argue that its dependence on the angle $ is weak, 
and a constant V c (d) = 0.26J was chosen in Ref. [98[ 
in order to satisfy the experimental constraint that the 
C-AF and G-AO order appears almost simultaneously 
in LaVOs [9l| . The experimental value T^f — 143 K 
for LaV03 [91( was fairly well reproduced in the present 
model taking J = 200 K. The functional dependence of 
the remaining two parameters {E z (i}), V a b(^)} on the tilt 
ing angle 1? was derived from the point charge model [9 
using the structural data for the -RVO3 series [lOll [lO 
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— one finds: 

E z {$) = Jv z sin 3 i? cos tf, (35) 
V ab (&) = J v a b sin 3 $ cos d . (36) 

Finally, the effective coupling to the lattice distortion 
has to increase faster with the increasing angle and 
the following dependence was shown [98| to give a sat- 
isfactory description of the phase diagram of the -RVO3 
perovskites: 

9cS (0) ~ J v g sin 5 d cos . (37) 

Altogether, magnetic and orbital correlations described 
by the spin-orbital model (|19[) . and the magnetic T/vi 
and orbital Too transition temperatures depend on three 
parameters : { v z , v & , v g } . 

Due to the proximity of both orbital and magnetic 
phase transitions in the i?V03 perovskites, it is crucial 
to design the MF approach in such a way that the spin- 
orbital coupling is described beyond the factorization of 
spin and orbital operators. On the one hand, the correct 
MF treatment of the orbital and magnetic phase transi- 
tions in the -RVO3 vanadates requires the coupling be- 
tween the on-site orbital, (t z )g = \\( T l — T j)\' an d sphr 
order parameters in the C-AF phase, (Sf)c, as well as a 
composite (S*t?) order parameter, similar to that used 
for the i?Mn03 perovskites [78j . On the other hand , the 
on-site MF theory including the above coupling |l04j does 
not suffice for the i?V03 compounds as the orbital singlet 
correlations (n-Tj) on the bonds (ij) along the c axis play 
so crucial role in stabilizing the C- AF p hase [79| and the 
orbital fluctuations are important |l05j . Therefore, the 
minimal physically acceptable approach to the present 
problem is a self-consistent calculation for a bond (ij) 
along the c axis, coupled by the MF terms to its neigh- 
bors along all three cubic axes (98|. This procedure, with 
properly selected model parameters, was shown to be suc- 
cessful in reproducing the experimental phase diagram of 
Ref. [92|]. One finds that indeed the orbital order occurs 
below a higher temperature than the magnetic one in 
the -RVO3 perovskites to the left from LaV03, i.e. with 
smaller ionic radius tr. 

As presented in Ref. |98[, the remarkable dependence 
of both spin T/vi and orbital Too transition tempera- 
ture in the -RVO3 perovskites follows from the respective 
changes in the orbital correlations with decreasing tr. 
First, the singlet correlations are drastically suppressed 
from LaV03 towards LUVO3. Second, the increase of 
orbital intersite interactions due to the JT term (131)1) . 
induces steady increase of the orbital temperature Too 
with decreasing tr. Finally, while (rf ) ~ 0.03 is rather 
weak in LaV03, it steadily increases along the i?V03 
perovskites when tr decreases, and finally it becomes 
as important as the orbital order parameter itself, i.e. 
(rf ) ~ (t*)q. Note that in the entire parameter range 
the latter order parameter is substantially reduced from 
the classical value (r z )G iinax = | by singlet orbital fluctu- 
ations in the entire parameter regime, being (rf )g ~ 0.32 
and 0.36 for LaV03 and LUVO3, respectively. 



It is quite remarkable that the above changes in the 
orbital state modify the magnetic exchange constants 
{Jab,Jc\ along both nonequivalent cubic directions, see 
Eq. (|24l) . and thus the value of T/vi is reduced with de- 
creasing tr. Note that the superexchange energy J does 
not change, so the entire effect stems from the orbital 
correlations |9Sj . This also implies that the width of the 
magnon band given at T = by Wc-af = 4(J a 6 + \J C \) 
is reduced by a factor close to 1.8 from LaV03 to YVO3, 
in agreement with surprisingly low magn on energies ob- 
served in the C-AF phase of YV0 3 [l06| . 

Summarizing, the microscopic model (|32[) describes 
gradual changes of the orbital and magnetic correlations 
under the coupling to the lattice which suppresses orbital 
fluctuations generated by virtual charge fluctuations re- 
sponsible for the spin-orbital superexchange. It provides 
an almost quantitative understanding of the systematic 
experimental trends for both orbital and magnetic tran- 
sitions in the HVO3 perovskites [98[ , and is able to repro- 
duce the observed non-monotonic variation of the orbital 
transition temperature Too for decreasing tr. However, 
the theoretical description of the magnetic transition to 
the G-AF phase at Tv 2 , which occurs for small tr [9lj |. 
remains to be addressed by future theory. More examples 
of spin-orbital entanglement in the field of the perovskite 
vanadates are shortly discussed in the next two subsec- 
tion. 



C. Optical spectral weights for LaVC>3 

As a second example of spin-orbital entanglement in 
the cubic vanadates at finite temperature we discuss 
briefly the evaluation of the optical spectral weights 
from the spin-orbital superexchange for LaV03, following 
Eq. (|23p . First we rewrite the superexchange operator 
H^(ij) for a bond (ij) \\ 7, contributing to operator T-Lj 
(fT9| . as a superposition of alfdj ^ g? 3 g?] charge excita- 
tions to different spin states in upper Hubbard subbands 
labelled by n [13], 

H^{ij)=Y,H { :} j . (38) 

n 

One finds the superexchange terms H„ \j for a bond (ij) 
along the c axis [8(| , 




and H n f ' for a bond in the ab plane, 
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n,ij 
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When the spectral weight is evaluated following Eq. 
it is reasonable to try first the MF approximation and to 
separate spin and orbital correlations from each other. 
The spectral weights require then knowledge of spin cor- 
relations along the c axis and within the ab planes (|28[). 
as well as the corresponding intersite correlations (t| • fj) 
and (t" t") with a — x,y, z. From the form of the above 
superexchange contributions one sees that high-spin exci- 
tations H^Jj support the FM coupling while the low-spin 

ones, and H-$}j, contribute with AF couplings. 

The low-energy optical spectral weight for the polar- 
ization along the c axis decreases by a factor close 
to two whe n the temperature increases from T ~ to 
T = 300 K |l07l j — this change is much larger than the 
one observed in LaMnOs (86|. However, the theory based 
on the MF decoupling of the spin and orbital degrees of 
freedom gives only a much smaller reduction of the weight 
close to 27%, and has no chance to explain the experi- 
ment as the maximal possible reduction of K ^ found for 
s c = in the limit of T — > oo amounts to 33% [28}. Note 
that both spin and orbital intersite correlations change 
in the temperature range < T < 300 K used in ex- 
periment, but this variation is clearly not sufficient to 
describe the experimental data. 

In contrast, when a cluster method is used to determine 
the optical spectral weight from the high-spin superex- 
change term (|39l) by including orbital as well as joint 
spin-and-orbital fluctuations along the c axis, the tem- 
perature dependence resulting from the theory follows 
the experimental data [80] . This may be considered as a 
remarkable success of the theory based on the spin-orbital 
superexchange model derived for the RVO3 perovskites, 
an the proof that spin-orbital entangled states contribute 
in a cruciat way in the finite temperature regime. In ad- 
dition, the theoretical calculation predicts that the low 
energy spectral weight is low along the c axis. The spec- 
tral weight in the ab planes behaves in the opposite way 
— it is small at low energy, and high (but not as high as 
the low-energy one for the c axis) at high energy. This 
weight distribution and its anisotropy between the c and 
ab directions reflects the nature of magnetic correlations, 
which are FM and AF in these two directions. A more 
precise comparison of these theory predictions for the ab 
polarization is not possible at present, but we expect that 
future experiments will also confirm them. 



sition at T N2 = 77 K from the G-AF to the C-AF spin 
order with rather exotic magnetic properties, found in 
YVO3 [108| . This magnetic transition is unusual and 
particularly surprising as the staggered moments are ap- 
proximately parallel to the c axis in the G-AF phase, and 
reorient above Xjv2 to the ab planes in the G-AF phase, 
with some small alternating G-AF component along the 
c axis. First, while the orientations of spins in G-AF and 
G-AF phase are consistent with the expected anisotropy 
due to spin-orbit coupling |l09j |. the observed magneti- 
zation reversal with the weak FM component remains 
puzzling. Second, it was a lso established by neutron 
scattering experiments 106] that the scale of magnetic 



excitations is considerably reduced for the G-AF phase 
(by a factor close to two) as compared with the exchange 
constants deduced from magnons measured in the G-AF 
phase. In addition, the magnetic order parameter in the 
G-AF phase of LaVC>3 is strongly reduced to ~ 1.3/is, 
which cannot be explained by rather small quantum fluc- 
tuations in the G-AF phase [99| ■ Finally, the G-AF phase 
of YVO3 is dimerized. Until now, only this l ast featur e 
found a satisfactory explanation in the theory [llOl lllll j , 
see below. 

We remark that the observed dimerization in the 
magnon dispersions may be seen as a signature of entan- 
glement in excited states which becomes active at finite 
temperature. The microscopic reason of the anisotropy 
in the exchange constants 



Jcl = Jc(l + 5 S ) , J C 2 = Jc(l - o~s) , 



(45) 



is the tendency of the orbital chain to dimeriz e, ac tivated 
by thermal fluctuations in the FM spin chain [lllj which 
support dimerized structure in the orbital sector. As a re- 
sult one finds alternating stronger and weaker FM bonds 
along the c axis (|4"5"|) in the dimerized G-AF phase (with 
S s > 0). The observed spin waves may be explained by 
the following effective spin Hamiltonian for this phase 
(assuming again that the spin and orbital operators may 
be disentangled which is strictly valid only at T = 0): 

H s = Jo {1 + (-!)**-}& -3+1 

(»,»+!> I|c 

+ Jab ]T S t ■ Sj +K Z J2 {S*f . (46) 

(ij)\\ab i 

Following the linear spin-wave theory the magnon dis- 
persion is given by 



1 



u±Qt) = 2\[2J ab + \J c \ + -K s ± JcVk 



with 



1/2 



(2Ja67k) 2 , 

(47) 



D. Peierls dimerization in YVO3 

The third and final example of the spin-orbital entan- 
glement at finite temperature in the family of vanadate 
perovskites is the remarkable first order magnetic tran- 



7k = 2 ( cosfc z + cos k y ) 
r]k = cos 2 k z + <5 2 sin 2 k z 



(48) 
(49) 



The single-ion anisotropy term oc K z is responsible for 
the gap in spin excitations. Two modes measured by 
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neutron scattering [106j are well reproduced by w±(k) 
obtained from Eq. P?]) using the experimental exchange 
interactions: J a b — 2.6 meV, J c = —3.1 meV, S s = 
0.35. We note that a somewhat different Hamiltonian 
with more involved interactions was introduced in ref. 
jl06l |. but the essential features seen in the experiment 
are well reproduced already by the present effective spin 
exchange model H s , see Eq. PB")) . 

The observed dimerization in the magnon spectra in 
YVO3 motivated the search for its mechanism within 
the spin-orbital superexchange model. Dimerization of 
AF spin chains coupled t o ph onons is well known and 
occurs in several systems Ill2j. The spin-Peierls transi- 
tion discovered in CuGeC>3 |l!3| led to renewed interest 
in the dimerization instability of the AF spin chains. In 
the spin-orbital model for the -RVO3 perovskites a sim- 
ilar instability might also occur without the coupling to 
the lattice when Hund's exchange is sufficiently small. 
In particular, the ground state at r\ — may be approx- 
imated by the dimerized chain with strong FM bonds 
alternating with the AF ones, if such chain s are coupled 
by AF interactions along the a and b axes [l!4{ (the ID 
chain would give the entangled disordered ground state 
as described in Sec. IV A|) . 

At realistic values of 77 > 0.10 the C-AF order with 
FM chains along the c axis is found in the ground state 
[79} . Numerical studies performed at finite temperature 
have shown that periodic dimerization of the magnetic 
exchange exists in a certain finite temperature range, 
while the ground sta te is the fully polarized and uniform 
FM state jl09l IllOj . These findings served as a moti- 
vation to investigate the mechanism of the spin-Peierls 
dimerization in FM spin chains. The microscopic ID 
model which stands for the situa tion encountered in the 
C-AF phase of YV0 3 reads [ml ]: 

h St = jJ2 (SvS i+ i + 1) (n-n+x + \-ih\ (50) 

where is stands for the contribution due to the high- 
spin states proportional to the Hund's exchange (|T8|) and 
stabilizes FM spin order. While the spin and orbital op- 
erators are disentangled in the FM ground state, one may 
consider a coupled FM spin chain to an orbital chain with 
interactions which favor the AO order, as realized in the 
C-AF phase. The exchange interactions along the spin 
(orbital) chain depend on the orbital (spin) correlations, 
and their modulation may be described by 6 S and S T pa- 
rameters (|45l) . They can be found from a self-consistent 
solution of the coupled MF equations for spin and or- 
bital correlations, and one finds indeed dimerized spin 
and orbital chains in a finite range of temperature [11 lj . 

Summarizing, spin-orbital entanglement in the excited 
states is also responsible for the exotic magnetic proper- 
ties of the C-AF phase of YVO3. They arise from the 
coupling between the spin and orbital operators which 
triggers the dimerization of the FM interactions as a man- 
ifestation of a universal instability of FM chains at finite 
T, which occurs either by the coupling to the lattice or 



to purely electronic degrees of freedom 111| . This latter 
mechanism could play a role in many transition metal 
oxides with (nearly) degenerate orbital states. 



VI. COEXISTING CHARGE AND ORBITAL 
ORDER 

The first step towards understanding the doped sys- 
tems with orbital degrees of freedom is the question con- 
cerning possible QP states deciding about coherent hole 
propagation in the orbitally ordered background. As 
discussed in Sec. [TT1 a single hole doped into the AF 
background as in Cu02 planes of La2Cu04 may propa- 
gate through the lattice because it couples to quantum 
spin fluctua tions and becomes dressed with a "cloud" 
of magnons jll5| . This results in the new energy scale 
oc J in place of the hole hopping t, so the hopping is 
strongly renormalized. The QP which forms after the 
hole is doped in the AF background is called a spin po- 
laron 36] . A more complex situation can occur in the sys- 
tems with partly filled degenerate orbitals, where a doped 
hole may not only couple to magnons but also couples to 
crystal- field excitations [l4| • In addition, QP states with 
higher spin states may occur, as for instance a triplet QP 
in case of an S = 1/2 antiferromagnet doped by a sin - 
gle electron in the orbitally degenerate background [l 16l j . 
This motivates two questions in the theory: (i) whether 
orbital excitations could couple as well to the moving 
hole and generate a new energy scale, as the magnons 
do, and (ii) whether spin-orbital entanglement has any 
important consequences for the hole dynamics. Both of 
them were addressed in the orbital t-J model for e g elec- 
trons 117], and in the analogous models for i 2g orbitals 
developed recently, see below. 

Two situations with a hole doped into an AO ordered 
background were considered in th e past: (i) a hole doped 
into an ab plane of LaMnOs [l 18j j which has an AO order 
of e g orbitals in the ground state, and (ii) a hole doped 
into an ab pla ne w ith an AO order of ti g orbitals and 
FM spin order [l!9j |. as realized for instance in Sr2V04. 
In the first case it was shown that the orbitons have 
in general a gap and have a lower dispersion than the 
magnons. Therefore, the quantum effects are weak but 
but a hole can move by interorbital hopping processes. 
While the constraint of creating no double occupancies 
has to be obeyed along the hole hopping, the bandwidth 
is strongly renormalized w ith r espect to that suggested 
by the LDA+C/ approach |l!8j |. Such interorbital hop- 
ping processes are absent in the t^ g ordered background 
with alternating yz/zx orbitals in an ab plane, and due 
to the specific t2 g orbital symmetries the o rbitons are dis- 
persionless. Thus the string picture 120] dominates the 
character of the ti g orbital polarons even more than in 
the case of systems with e g orbital degrees of freedom. 

An intriguing question in this context addressed only 
recently is whether spin quantum fluctuations can still 
contribute to the QP spectral properties when both types 
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of order, spin and orbital, alternate in an ab plane, as for 
instance in the C-AF phase of -RVO3 perovskites. A cru- 
cial observation for the spectral properties of a hole doped 
into the entangled AF / AO background is a simultaneous 
excitation of a magnon and an orbiton when a hole moves 
by a single step in the lattice |l2l| . This dominates the 
behavior of the hole doped in such an entangled state, 
because the orbitals confine the hole motion by forcing 
the hole to retrace its path which implies that the hole 
motion by its coupling to the quantum spin fluctuations 
is prohibited. Thus, the string-like potential which acts 
on the hole is induced by the orbitals although it has a 
joint spin-orbital character. Hence, this important fea- 
ture of the orbitally induced string formation could be 
understood as a topological effect. This happens even 
if the energy of the orbital excitations is turned to zero, 
i.e., when the hole moves in the orbital sector incoher- 
ently. Hence, the mere presence of orbitals is sufficient 
to obtain the (almost) classical behavior of a hole doped 
into the ground state with AF/AO order. This result, 
in connection with the fact that the mother-compound 
of the superconducting iron-pnictides shows a variety of 
spin-orbital phenomena [l22j | , suggests that further inves- 
tigation of the hole propagation in spin-orbital systems 
is a fascinating subject for future studies. 

The properties of doped fl1_.rSr.rVO3 systems are 
puzzling and it is not understood until now why (i) 
Lai-^Sr^VOs is insulating in a broad range of doping 
below x c = 0.18, and (ii) why the AF order survives 
even for x > x c when the system beco mes metallic and 
looses the AO order in the ab planes |l23f . The grad- 
ual changes of the optical conductivity under increasing 
doping demonstrate that the anisotropy between the ab 
and c direction decreases, but surprisingly is not com- 
pletely lost even in the metallic regime. The differences 
observed in the optical conductivity and raman scatter- 
ing spectra between Lai-^Sr^VOs and Yi-^Sr^VOa sug- 
gest that the orthorhombic lattice distortion plays also 
here a very important role and influences the hole dy- 
namics |l24j . It has been argued that the C-AF phase 
is more robust under hole doping [l 251 ] and survives in 
a broad doping range [l26j . but a complete understand- 
ing of doped vanadate perovskites awaits a more careful 
theoretical study. 

Doped i?i_ :E (Sr,Ca) 2 :Mn03 systems are studied much 
longer and they are better understood. The FM metallic 
state is induced by doping via the double exchange mech- 
anism U. I wh ich was also formulated for degenerate e g 
orbitals 127j , and the phase diagrams of the doped per- 
ovskite syst ems show a remarkable sequence of magnetic 
phases [l28j . from the A-AF phase, through the insulat- 
ing and metallic FM phase, towards the C-AF phase and 
G-AF phase in the highly doped regime. Similar (but 
not the same) sequence of magnetic phases was r eported 
in the bilayer La2-2xSri_|_2xMti207 systems [l29| . As in 
the -Ri^Sr^VOa perovskites, also in the layered systems 
the orbital ordered (or liquid) state determines whether 
the intersite spin correlations are AF or FM, as shown 



for the monolayer |130j and for the bilayer [131| using a 
t- J-like model which includes orbital degeneracy. These 
model calculations illustrate as well the complementarity 
of spin and orbital order expressed by the Goodenough- 
Kanamori rules [72j. In these systems the short-range 
charge order gradually develops with increasing doping 
in the realistic parameter regime [132| . However, more 
complete models including the charge transfer physics are 
necessary to describe the features observed in the optical 
spectra, as for instance in insulating LaSrMn04 133]. 

Although there is no complete understanding of the 
phase diagram and in particular of the mechanism of 
the metal-insulator transition which leads to the colos- 
sal magnetoresistance until now, a lot of progress could 
be made using model Hamiltonians. It was recognized 
that the orbital dege nerac y plays a crucial role both in 
the double exchange 127j and for the interactions with 
the lattice due to the JT effect jl34| , and phase diagrams 
which resemble the qualit ative behavior of the doped 
manganites were obtained fl35]. However, a more real- 
istic treatment requires also electron correlations among 



e„ electrons which are more difficult to implement 
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contrast to the nondegenerate Hubbard model, the or- 
bital Hubbard model for FM manganites does not show 
an instability towards the orbitally polarized FO state 
and o ne fin ds instead the disordered orbital liquid ground 
state [136]. This concept was crucial in explaining the 
doping dependence of the st iffne ss constant in the FM 
Lai_2a;Sr a ;Mn03 manganites 137J, but for a quantitative 
explanation both the double exchange due to correlated 
e g electrons and the superexchange due to ti g core spins 
had to be included |l38l|. This approach had also a re- 
markable success [139| in explaining the observed magnon 
dispersion and the doping dependence of the magnetic 
exchange constants in the La2-2a;Sri_|_2 a: Mn207 systems, 
including the observ ed p hase transition from the FM to 
the A-AF structure [l40j . 

As expected, the orbital order melts in general when 
the manganites are doped, there are cases when a dif- 
ferent type of orbital order coexisting with charge or- 
der emerges again at half doping. The famous case is 
the so-called (charge exchange) CE phase in half-doped 
(x = 0.5) manganites [7l|, where the two-sublattice 
charge order coexists with orbital order on the sites with 
the majority of e g electron charge, and the FM zig-zag 
chains staggered in ab planes. Although the double ex- 
change provides some argume nts j ustifying the stability 
of this complex type of order [l4lj which competes with 
the FM phase in the relevant parameter regime and wins 
for sufficiently large and realistic AF su perexchange be- 
tween the S = 3/2 t2 g core spins [l42j . the problem is 
subtle and the range of parameters with the CE phase in 
the ground state is rather narrow. In particular, this 
phas e is destabilized by intersite Coulomb interaction 
[13l| , and the JT distortions play an important role in 
stabilyzing it [l43j . We emphasize again that AF inter- 
actions between ti g electrons are small — for this case 
the CE phase was found in the charge-ordered phase us- 
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ing a finite-temperature diagonalization technique [l44j . 
The mechanism invoked there to stabilize the CE phase 
is subtle and employes the cooperative JT interaction 
between next-nearest Mn 3+ neighbors mediated by the 
breathing mode distortion of Mn 4+ octahedra and dis- 
placements of Mn 4+ ions. It is worth noting tha t the 
topological phase factor in the Mn-M n ho pping [l45j | 
leading to gap formation in ID models [l4l| for the CE 
phase, as well as the nearest neighbor JT coupling, are 
not able to produce t he o bserved zigzag FM chains for the 
realsitic parameters [l44j . Recent x-ray structural analy- 
sis of Pr .5Cao.5Mn03 and Euo.sCai.sMnC^ suggest that 
the orbital shape and the charge disproportionation are 
sensitive to the dimension of Mn-0 network [l46j | , which 
together with the possible different role of the JT effect in 
both compounds poses new interesting questions in the 
theory. 

The controversy about the nature of the charge order 
in this phase which arose due to signa tures of Zener po- 
larons observed in the neutron data [1471 ] seems to be 
resolved now in f avor of the more conventional picture 
of zig-zag chains |148j . Melting of this composite order 
with increasing temperature is fascinating and the mag- 
netic order disappears first, giving coexisting charge and 
orbital order in the intermediate temperature regim e, be- 
fore both melt resulting in a disordered phase [149]. This 
suggests that the link between (weak) charge and orbital 
order is particularly strong here, similar as in the mag- 
netite below the Verwey transition [TTj] . Recently, charge 
ordered AF phase was also reported in Lai.sSro.sCoOa, 
and preliminary theoretical concepts in t he fr amework of 
spin-orbital physics were also presented [l50j |. 

Stripe phases appear also in the doped systems with 
active orbital degrees of freedom, but are qualitatively 
different from the ones observed in the cuprates . see Sec. 
III1 Stripe order was found in doped manganites [15lj and 
also discovered in doped L^-^Sr ^NiC ^ nickelates about 
the same time as in the cuprates [152 ]. However, in con- 
trast to the cuprates the stripes in La 5 / 3 Sr 1 / 3 Ni04 are 
dia gona l and contain one (and not half) hole per unit 
cell [153J . Intriguing features seen in the spin excitation 
spec tra of La2Ni04+5 nickelates were reported recently 
[l54l j which suggest that the inward dispersion, seen also 
in cuprates, has a common origin in stripe phases. Simu- 
lations performed within the LDA+t/ approach suggest 
that a subtle interplay between the charge and spin order 
and octahedral distor tions is essential for the formation 
of an insulating state [155| . Preliminary HF calculations 
emphasize the importance of orbital degeneracy i n th e 
case of nickelates for the realistic e g hopping model (l56j j , 
where one finds indeed that diagonal stripes with the ob- 
served filling (of one hole per unit cell) are more stable 
than other phases, in contrast to the predictions of the 
degenerate Hubbard model with diagonal hopping (which 
does not agree with experiment). It remains a challenge 
for the theory to develop a more complete theory of the 
stripe phases in the nickelates, including the electron cor- 
relations and the coupling of e g electrons to the lattice 



distortions, and to understand better differences between 
the stripe phase in the nickelates and in the cuprates . 



VII. SUMMARY AND OPEN PROBLEMS 

Summarizing, charge order is common in doped tran- 
sition metal oxides and arises even in absence of intersite 
Coulomb interactions while the kinetic energy of doped 
charges (holes or electron) competes with the magnetic 
superexchange. Under such circumstances stripe phases 
realized in the cuprates, nickelates, and manganites are 
favored as then the two above energies are optimized si- 
multaneously in the domain walls and in the magnetic 
domains between them. When also orbital degrees of 
freedom are active, the charge order is accompanied by 
certain, usually weak, orbital order. Good examples of 
this composite type of charge-and-orbital order are the 
low temperature phase of Fe3 04, and the CE phase in 
the half-doped manganites. 

In the perovskite lattice the orbital order is stabilized 
easier in the correlated insulators with orbital e g degrees 
of freedom, as then the superexchange and the interac- 
tion with the lattice act supporting each other [78| , than 
in the ones with the ti g active orbitals, where the inter- 
actions with the lattice in general compete with the su- 
perexchange (98|. As a result, the orbital and magnetic 
transition occur independently from each other and at 
quite different temperatures in the e g systems, while the 
case of the .RVO3 perovskites is an example of the prox- 
imity and interplay of the magnetic and orbital phase 
transition. Both these different situations were success- 
fully described within spin-orbital superexchange models 
with added interactions with the lattice. 

A qualitative difference between the orbital order in e g 
and ti g systems is that the orbital order is more robust 
in e g case and may be usually treated by classical (MF) 
approaches, while t2 g orbitals may easier fluctuate and 
thus couple also easier to the spin degrees of freedom. It 
is for this reason that composite spin-orbital fluctuations 
occur in correlated titanates and vanadates insulators. 
Although such fluctuations are quenched in the ground 
state of these systems for realistic parameters, they de- 
velop at increasing temperature due to the presence of 
excited stated with spin-orbital entanglement. For in- 
stance, such composite spin-orbital fluctuations are re- 
sponsible for the temperature dependence of the optical 
spectral weights in LaV03 [8(| and trigger spin-orbital 
dimerization in the C-AF p hase of YVO3 in the inter- 
mediate temperature regime [l09| ]. Similarly interesting 
orbital ordered states are also found in the perovskite 
ruthenates — as an example we mention here the puz- 
zling low temperature electronic and structural behavior 
recently discovered in PbRu03 [l57j ] . 

The microscopic mechanism of melting of the orbital 
order in doped systems is a very challenging problem 
in the theory and could not be understood until now. 
The main difficulties follow from disorder and the ne- 
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cessity of using Monte Carlo techniques. In this way it 
could be concluded that a strong competition between 
the FM metallic and t he AF charge-ordered insulating 
states takes place [l58| |. This competition influences the 
transport properties and leads to short-range spin and 
charge correlations which evolve with time. It is intrigu- 
ing to what extent this complex situation in the metallic 
phase influences the magnetic excitations of the system. 
The magnetic excitations in the metallic FM phase of sev- 
eral doped manganites soften at the zone boundary and 
may be described by the Heisenberg model with the near- 
est neighbor Ji and fourth-nearest neighbor J4 exchange 
coupling |l59j | . Its microscopic origin is controversial and 
two different concepts were proposed to explain the ex- 
perimental observations: (i) quantum fluctuations of the 
slanar x 2 — y 2 orbitals associated with the A-AF phase 
163, and (ii) the (3z 2 — r 2 )-type orbital fluctuations 
159]. It was also shown that the ratio of J4/J1 changes 
along the Ai_ x A' x Mn03 manganites (with A and A' be- 
ing the rare-earth and alkaline rare-earth ions), while the 
stiffness constant is almost universal and has only very 
weak dep ende nce on the chemical composition for a fixed 
doping x [l6l| . Both theoretical models have difficulties 
to explain the experimental data — an incorrect disper- 
sion along the (111) direction follows from the first one 
[I6dj . while the second one predicts a spectacular dop- 
ing dependence of J4/J1 which is not observed. Thus, 
in spi te of its remarkable success in the simplest situ- 
ation 138], the complete theoretical explanation of the 
magnon dispersion in the metallic FM manganites within 
the orbital liquid state remains a challenging problem 
in the theory. Even more puzzling are the magnons in 
the insulating FM phase, where several bra nche s with a 
staircase-like spectrum have been observed |162j |. 

Recent studies of the phase diagram of the i?MnC>3 
manganites include the effect of orthorhombic distortions 
within theoretical models with spin superexchange and 
the Dzyaloshinsky-Moriya interaction responsible for the 
multiferroic behavior [1631 ) . The usual approach so far 
is to develop an effective spin model including the spin - 
lattice coupling, leading to the electric polarization |l64j . 
It seems that an explicit treatment of the orbital degrees 
of freedom could provide a better understanding of the 
observed phenomena. 

Another active direction of research in the field of 
transition metal oxides is the search for novel quantum 
phenomena, including more examples of quantum spin- 
orbital entanglement. They could be found in frustrated 
lattices, and here we mention briefly only the triangu- 
lar lattice. An interesting case and good candidate for 



a spin-orbital liquid might be LiNi02 with a triangular 
lattice of Ni 3+ ions (d 7 configuration with 5=1/2 spins) 
and no magnetic or orbital order down to very low tem- 
peratures. First, it was a rgued that a model based on 
symmetry arguments jl65| ] , characterized by a large num- 
ber of low-lying singlets associated to dimer coverings of 
the triangular lattice, could explain the absence of any 
type of ordered phase in L1N1O2. Second, it was shown 
that excited states on oxygens along the 90° bonds are 
crucial in the superexchange and they change the balance 
between different terms in the Hamiltonian, making the 
orbital interactions stronger than the spin ones [166J. In 
any case, interplane JT coupling seems to be too weak 
in L1N1O2 to stabilize the orbital long-range order, and 
the microscopic reason of disorder could be alone due to 
strongly frustrated orbital interactions on the triangular 
lattice in (111) planes [l66| . which resemble the compass 
model. A completely different situation is encountered 
in the d 1 spin-orbital model with active three ti g orbitals 
on the triangular lattice, as realized in NaTi02, where a 
spin-orbital disordered liquid state is more likely [l67j . 

As a final remark, we would like to mentio n rec ent ex- 
perimental studies of Ni-based superlattices 168} . They 
stimulated progress in the theory which predicts that, 
in analogy to t he g ain of kinetic energy in the lay- 
ered manganites [1691 ]. the correlated e g electrons in the 
MO2 planes develop a planar {x 2 — y 2 )-like orbital or- 
der in LaNi0 3 /LaM0 3 superlattices (with M=A1, Gd, 
Ti) |l70| . It may be expected that future studies of the 
systems of reduced dimensionality will provide more un- 
expected properties in the near future, and could lead to 
developing functional materials, using both charge and 
orbital degrees of freedom. 
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